arXiv:1501.02956v2 [astro-ph.SR] 5 May 2015 


SIMULATIONS OF THE SPATIAL AND TEMPORAL INVARIANCE IN 
THE SPECTRA OF GRADUAL SOLAR ENERGETIC PARTICLE 

EVENTS 


Yang Wang^ and Gang Qin^ 

ywang@spaceweather.ac.cn; gqinOspaceweather.ac.cn 
Received_; accepted_ 


'State Key Laboratory of Space Weather, Center for Space Science and Applied Research, 


Chinese Academy of Sciences, Beijing 100190, China 





ABSTRACT 


The spatial and temporal invarianee in the speetra of energetie partieles in the grad¬ 
ual solar events is reprodueed in the simulations. Based on a numerieal solution of the 
foeused transport equation, we obtain the intensity time profiles of solar energetie 
partieles (SEPs) aeeelerated by an interplanetary shock in the three-dimensional inter¬ 
planetary space. The shock is treated as a moving source of energetic particles with a 
distribution function. The time profiles of particle flux with different energies are cal¬ 
culated in the ecliptic at 1 AU. According to our model, we find that shock acceleration 
strength, parallel diffusion and adiabatic cooling are the main factors in forming the 
spatial invariance in SEP spectra, and perpendicular diffusion is a secondary factor. In 
addition, the temporal invariance in SEP spectra is mainly due to the effect of adiabatic 
cooling. Eurthermore, a spectra invariant region, which agrees with observations but 
is different than the one suggested by Reames and co-workers, is proposed based on 
our simulations. 

Subject headings: Sun: activity — Sun: coronal mass ejections (CMEs) — Sun: particle 
emission 
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1. INTRODUCTION 


Solar energetic particl e (SEP) events can r oughly be divided into two categories: impulsive 


events and gradual events (iReames 


19951 


19991). The impulsive events, with the characteristics 


of low intensity and short duration, are produced by solar flares. Gradual events, usually lasting 
longer and having h igh intensity, are r elated to the shocks driven by interplanetary coronal mass 

(120061) investigated the radial and longitudinal dependence of 


ejections (ICMEs). 


Lario et al. 


4-13 and 27 - 37 MeV proton peak intensities and fluences measured within 1 AU. They found the 
peak intensities and fluences of SEP events can be approximated by j = jQr~" exp 
where j is either the peak intensity or the fluence, r is the radial distance of the spacecraft, 0 is the 
longitudinal angular distance between the footpoint of the observer’s field line and the region of 
the SEP source, and 0o is the centroid of the distributions. Eurthermore, the radial dependence of 
peak inte nsities and fluence s of SEP events have been simulated with a focused-diffusion transport 


equation (|Lario et alJ 


2007h . 


shocks: some aut 

lors ( 

Heras et a 

.1992, 

1995 

: Kal 

enrode & Wibberenz 

1997; 

Kallenrode 

2001; 

Ng et al. 

2003; 

Wang et al. ^ 

>012; 

Oin et al. 

2 OI 3 I) adopted a ‘ 


)V CME driven 


Eario et al. 


1998 


to treat the shock as a moving source, and SEPs are injected at the shock with an assumed 


injection strength, while a few other studies include the acceleration of SEPs by shocks 


1983 


Gordon et al 


1999 


Zank et al 


20001 : 


Ei et al. 


2003 


Rice et al. 


20031: 


Sokolov et al. 


Lee 


2004 


Li et al. 


2005 


Kota et al. 


200.51: 


Tylka & Lee 


2006 


Zuo et al 


2011 


Zuo et al. 


2013). In these 


models, three important effects of acceleration and propagation mec hanisms have been involved. 


The first effect is the acceleration process by the CME-driven shock. 


Zank et al. 


(120001) modeled 


the evolution of a CME-driven shock based on an “onion shell” mode]j_and this model has been 


furtherly developed in a number of papers (iLi et al. 


2003 


Rice et al. 


2003 


Li et al. 


2005b. They 


used a magnetohydrodynamics (MHD) code to describe the evolution of the CME-driven shock 



































































































































-4- 


in the interplanetary spaee, and wave exeitation by streaming energetie partieles produeed at 


shoek is ineluded. Based on the model, the simulation result s ean sueeessfully exp 


fluxes and speetra in some multi-spaeeeraft observed events (IVerkhoglyadova et al. 


ain the SEP 


200 ^ 


he seeond effect is energetic particles to interact with Alfven wave self-consistently. 


2010 ). 


(2003 


Ng et al 


2012|) presented a model of particle transport including streaming proton-generated Alfven 


waves, and the amplification of the Alfven waves is determined by the anisotropy of particles. 
The particle diffusion coefficients can be calculated from wave intensity and wave growth rates. 
Their simulation results show a good agreement with the observed spectral slope a nd abundance 


ratios of heavy ions. T h e third effect is the rea listic geometry of CME and its shock dSokolov et al. 


2004; 


Kota et al. 


20051). 


Sokolov et al.l (l2004l) modeled particle acceleration and transport as CME 


driven shock wave propagating from 4 to 30 solar radius from the Sun. The realistic structures 
of CME and its shock are derived from a numerical solution of a fully three-dimensional MHD 
model. Their simulation results demonstrate that the diffusive shock acceleration theory can 
account for the increase of hundreds of MeV protons during the early stages of CME driven shock. 

Solar energetic particle events measured by multi-spacecraft help us to understand the 
processes of particle acceleration and transport in the heliosphere. In some gradual events, 
the SEP fluxes measured by widely separated spacecraft presen t similar intensities within a 


smal 


2010 


2 - 3 factor in different lat i tudes, longitudes or radius (IReames et a 


2013 


McKibben et al 


200 Ibl: 


Maclennan et al. 


This phenomenon w as firstly proposed by 


2001 


Eario et al 


2003 


1997 


; IReames 


Tan et al 


2009). 


Roelof et al. 


McKibben et al 


McKibbenI (Il972h . and wa s named “r e servoi r” by 


(Il992t). In order to interpret the reservoir phenomenon. 


gradients of flux, while 


McKibbenI (Il972h and 


(l2001b) involved an effect ive perpendicular diffusion to reduce the spatial 


Roelof et al. 


(Il992h suggested a diffusion barrier produced by ICMEs 


or shocks. The magnitude of magnetic field increases at the outer boundary of reservoirs, so 
that SEPs could be contained in the reservoirs for a long time. Eurthermore, the i nterplanetary 


magnetic field (IME) has been disturbed by ICMEs, SEPs could be redistributed. 


Reames et al. 
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(119961) shows that, in some gradual SEP even ts, the spectra are inva riant both in space an d time 


This discovery extended the original work of 


McKibbenI (119721) . In 


Reames et al. 


(119961) . they 


considered an expanding magnetic bottle of quasi-trapped particles between an ICME driven 
shock and the Sun. As the magnetic bottle expanding, the SEP fluxes gradually decrease as a 
result of parallel diffusion and adiabatic cooling. In this sense, the magnetic bottle plays a pivotal 
role in the decay phase of SEP event. 

In principle, the disturbances in the magnetic field caused by ICMEs can help the particles 
redistribute in space. However, the reservoir phenomenon cannot be simply explain ed as a result 


of the disturbances o 


( 19961) and [Reames 


IMP caused by ICMEs. Pirstly, in the redistribution process in 


Reames et al. 


(12013h . no explicit transport mechanism can reduce latitudinal, longitudinal. 


and radial gradients of SEP fluxes besides perpendicular diffusion. Secondly, in some SEP events, 
CMEs are not directly observed by the spacecraft, the reservoir phenomenon is also observed 


McKibben et al. 


2001al) . Thirdly, in 


ReamesI (Il999t) . when the observer is located at the eastern 


of the shock, the onset time of temporal invariance in the SEP spectra is earlier than the shock 
arrival time. These results are not consistent with that of an expanding magnetic bottle. 

The effect of perpendicular diffusion is important in the SEP fluxes especially when the 
observer is disconnected from the shock by IMP. During the time period March 1, 1979 to March 
11, 1979, a gradual SEP event has been detected by Helios 1, Helios 2, and IMP 8. The three 


spacecraft are located in the ecliptic near 1 AU, but at different longitudes. In the c 


of thi s SEP event, the reservoir phenomenon appeared (IReames et al. 


1997 




hase 


201 


2013h . In this event, the in-situ observation shows that an ICME was detected by Helios 1, but not 


by Helios 2 an d IMP 8. And the interplanetary shock was only observed by Helios 1 and 2, but 


not by IMP 8 (lEario et al. 


20061 : 


Reames 


2010h . According to the location of the three spacecraft. 


if the ICME was located behind the center of the shock front, then Helios 1 was located near to the 
center of shock, and Helios 2 and IMP 8 were located at the West frank of the shock. However, 
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the onset time of SEP fluxes observed by three spacecraft was very close. How could SEPs be 
detected by IMP 8 before it was connected to the shock by IMP? One possible answer is the effect 
of perpendicular diffusion which also possibly works in forming the reservoir phenomenon. 

However, perpendicular diffusion has been always a difficult problem for several decades. 
Observation results show various levels of perpendicular diffusion coeffici ents for different S EP 


events. Por example, ‘dropout’ phenomenon in the impulsive SEP event (iMazur et al. 


200011 


usually show reduced perpendicular diffusion, in order to reproduce the ‘dropout’ phenomenon 
in simulations, perp endicular diffusion coefficient k i should be seyeral order magnitude smaller 


than 


parallel one xn (IGiacalone et al. 


20001 : 


Guo & Giacalone 


2014 


Droge et al. 


20101 : 


Wang et al 


2014 . On the other hand, for some events, observation result s show that the perpendicular 


diffusion coefficients cou 


2003 


Dresing et al 


d be comparable to the parallel ones (IDwyer et alJ 


1997 


Zhang et al. 


2012|1 . In order to understand diffusion, many efforts were made the oretically . 

3 


^ ly ass uming energetic particles’ perpendicular and parallel diffusion do not interaction. 


Jokipiil 


1966t) developed the quasi-linear theory (QET). According to QET, perpendicular diffusion 


coefficient is usually much smaller than the parallel one. However, it is found that intera ction 


between par allel and perpend icular diffusion is important in theory (IKota & Jokipi 


simu 


ations (lOin et al. 


2002b 


at), so the non-linear guiding center (NEGC) theory (iMatthaeus et al. 


2000l) and in 


2003(1 is developed to describe perpendicular diffusion with the influence of parallel diffusion. 


which agrees with simulations much bet ter than QET. In addi tion, simulations show different 


levels of perpendicular diffusion, e.g., in 


Pin & Shalchil (120121) . the magnitude of /Cx/zcy could be 


as large as 10 ^ in some conditions, and as small as 10 in other conditions. 


Recently, 


Pin et al. 


(120131) proposed that the shock acceleration strength makes important 


contributions to the reservoir phenomenon, particularly in low-energy SEPs. In their simulations, 
the reservoir phenomenon is reproduced under a variety of condition s of shock accele ration 


strength and perpendicular diffusion. In this paper, as a continuation of 


Pin et al. 


(|2013[) . we study 
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the property of SEP spectra in the decay phase. We compute the time profiles of SEP flux which 
are accelerated by interplanetary shock. In section 2 we describe the SEP transport model and the 
shock model. In Section 3 we show the simulation results. In Section 4 we summary our results. 


2. MODEL 


2006 


n this work, we model the transport o 


Zhang et al. 


Pin et al. 


2013 


2009- 


)r6ge et al 


Zuo et al 


equation is written as dSki11inp 


2013 


12010 : 


SEPs following previous research(e.g 


de et al. 


Wang et al. 


1971 


Schlickeisei 


2011 


Zuo et al 


mA 


Pin et a' 


Wang et al. 


2012 


201^^. A three-dimensional focused transpo rt 


20021 : 


: Pin et al. 

o 

o 

<N 

Zhang et al. 

o 

o 


df / ^ 

^ = V-(K^.Vf)-[viub+V^ 


■’ dii\ 


+P 


\-l? 


A A 


V V^^-bb : yV^] + iu^bb : VF 


A A 


1 — Ip [ V I AA 

--+//^v-r'"-366: vr 


d± 

dp 

dp' 


( 1 ) 


where f(x,p,p,t) is the gyrophase-averaged distribution function; jc is the position in a 
non-rotating heliographic coordinate system; t is the time; p, p, and v are the particle pitch-angle 

A 

cosine, momentum, and speed, respectively, in the solar wind frame; is a unit vector along the 
local magnetic field; F™ = F^*^ r is the solar wind velocity; and L is the magnetic focusing length 
given by L = with Bq being the magnitude of the background magnetic field. The 

IMP is set as the Parker field model, and the solar wind speed is 400 km/s. This equation includes 
many important particle transport effects such as particle streaming along the field line, adiabatic 
cooling in the expanding solar wind, magnetic focusing in the diverging IMP, and the diffusion 
coefficients parallel and perpendicular to the IMP. 


1986 


Pin et al. 


The pitch angle diffusion coefficient model is set as (IBeeck & Wibberenz 
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2005h 


D, 




Dqvp^ ^ +^)(i 


( 2 ) 


where the constant Dq controls the magnetic field fluctuation level. The constant q is chosen as 
5/3 for a Kolmogorov spectrum type of the power spectra density of magnetic field turbulence 
in the inertial range. Furthermore, h = 0.01 is chosen for the non-linear effect of pitch-angle 


diffusion at // = 0 in the solar wind (|Oin & Shalchil 


20091 


2014). 


_ The parallel mean free path (MFP) Tn can be written as (IJokipiill 19661:iFfasselmannI1 1 9681: Earl 

I974I1 


fill 


3v 


— 


2\2 


D 


■dp. 


(3) 




and the parallel diffusion coefficient /cn can be written as /cn = vdy/S. 

The relation of the particle momen tum and the perpendicular diffusion coefficient is set as 


Potgieter & Moraal 

1985; 

Zhang 

1999) 




-)(l-bb 


(4) 


where Be is magnetic field strength at the Earth, B is the magnetic field strength at the location of 
particle, p is particle momentum, and a is set to 1/3. Different perpendicular diffusion coefficients 


could be obtained by altering kq. Note that we use t 
the parameters, e.g., a could be set as other values 1 


lis ad-hoc m odel for the purpose of simplicity. 


Zhang 


19991) . However, the variation of these 


parameters would not qualitatively change the results in this paper. There are some more complete 
models that are developed to des cribe the particle diffusion in magnetic tur b 


nonli near guiding center theory (IMatthaeus et al 


2014 . 


2003 


Shalchi et al. 


2004. 


)ulence , such as th e 


2010 ; 


Pin & Zhang 


We use a time-backward Markov stoc hastic process met hod to solve the tr ansport equation 


([U). The detail of method can be found in 


Zhang 

(1999 

) and 

Pin et al. 

(2006) 


injection on the shock is specified by boundary values. The boundary condition is chosen as 
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2013h 


follow ing form (iKallenrode & Wibberenz 


1997; 

Kallenrode 

2001; 

Wane et al. 

2012; 

Pin et al. 


Mr,e, (p,p,t) 
S{r,e, (p,p) 


a-Sir- v,t) ■ Sir,e,(p,p)- p ^ ■ ^(6, (p) 

(_) . exp ——— 

rc (pc{p) 


I 1 if \(p{e, ip)\ < (f), 

0 otherwise, 


(5) 


where the partiele are injeeted at r = v^t, Vs is shoek speed, vj = ro + n- Ar, with r = 0,1,2 • • • no- 
Ar is spaee interval between two ‘fresh’ injeetions, ro = 0.05 AU is inner boundary, is set to 
0.05. r is distanee between sun and shoek. (p is the angle between the eenter of shoek and the 
point at the shoek front where the partieles injeeted. The shoek acceleration strength is set as S for 
specifying the particles ejection. It changes with a power law in radial distance and exponential 
towards the flank of shock. ^ determines the spatial scale of shock front, cps is the half width of the 
shock. 


3. RESULTS 


The parameters used are listed in table [B unless otherwise stated. Note that the IMF is set 
as Parker spiral, and the disturbances of IMF behind the shock are ignored. The particle energy 
channels are chosen as 5 MeV, 10 MeV, 20 MeV, 40 Me V and 80 MeV. T he parallel mean free 


path depends on the momentum Tn ~ p^^^. According to 


Pin et al. 


(12013h . the Kj_/ki\ is set as 0.1 


in the ecliptic at 1 AU. Because the shock acceleration efficient decreases as the particle energy 
increasing, the acceleration strength parameters also change with the momentum: a ~ 

(pc ~ The observers are located in the ecliptic at 1 AU. 
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3.1. Temporal Invariance in the Spectra 

In Figure [U and [2l we plotted the fluxes of different energy channels in the cases with and 
without adiabatic cooling. In order to check the temporal properties of SEP spectra in the decay 
phase, we normalize different energy fluxes, so that the fluxes have similar values soon after all of 
them have reached peaks. The panels show the normalized fluxes observed in the ecliptic at 1 AU, 
but at different longitudes E60, E20, W20 and W60. The notations E60, E20, W20, W60 are short 
for East 60°, East 20°, West 20° and West 60°, respectively. East/West means the location of the 
observer is east/west relative to the center of shock. The vertical lines indicate the shocks passage 
of the observers. 

In Eigure[Il the adiabatic cooling effect is included in SEP propagation process. In the 
decay phase of SEP events, shock acceleration strength, adiabatic cooling, parallel diffusion, and 
perpendicular diffusion are the major factors to influence the flux behavior. In the four panels, the 
fluxes at all energies follow a similar trend, then the fluxes scatter slowly as time goes by. This is 
called temporal invariance in the spectra of gradual SEP event. In the E60 event, the fluxes at all 
energies start to follow a similar trend about one day before the shock passage of 1 AU. In other 
words, the onset time of the temporal invariance is earlier than the time of the shock passage of 
the observer. In the E20 and W20 events, however, the onset time of temporal invariance is close 
to the shock passage of the observers. In the W60 event, furthermore, the temporal invariance 
starts the latest, and actually it starts two days latter than the shock arrival. 

In Pigure|2l the adiabatic cooling is not included in the SEP propagation process. Without 
adiabatic cooling, shock acceleration strength, parallel and perpendicular diffusion are the major 
factors in the decay phase. Due to the different diffusion coefficients and shock acceleration 
strength for different energy particles, the fluxes decay with different ratios consequently. With 
higher energies, the fluxes decay much faster. In these cases, the temporal invariance does not 
exist in the decay phase. 
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Comparing Figure [U and Figure [2l the fluxes deerease mueh faster with adiabatic cooling. 
Because of adiabatic energy loss, particles have less energy when they are observed than that 
when they are released in the sources. In addition, since the source spectrum index is negative, 
the fluxes are lower with higher energies, so that the adiabatic cooling effect makes the SEP flux 
decreasing as time passes by. To sum up, the temporal invariance in the spectra results from the 
adiabatic cooling effect. 


3.2. Spatial Invariance in the Spectra 

In Figured the SEP fluxes are shown for three observers located at different longitudes, E20, 
W20, and W60. The upper panel shows the 5 MeV proton fluxes observed by the observers. We 
set two typical time intervals, interval A from 1.3 days to 1.5 days in rising phase and interval 
B from 6.9 days to 7.1 days in decay phase. In order to study the spatial variance in different 
phases, in lower left and right panels of PigurelH we plot the energy spectra observed in different 
longitudes in interval A and interval B, respectively. During the interval A, the spectra are different 
among the three observers. However, during the interval B, spectra are almost the same among 
This phenomenon, which is named spatial invariance in the spectra by 
I, results from the reservoir effect in different energy channels. 

Because the shock is a moving source in the interplanetary space, the peak intensity of SEP 
flux is mainly determined by the shock acceleration strength and parallel MEP In upper panel of 
Eiguredl at the peak time of flux for W20 (W60), the flux for W20 (W60) is close to that for E20. 
Eurthermore, the SEP fluxes decay in a similar ratio because of the effect of adiabatic cooling. 

At the same time, the latitudinal gradient in the SEP fluxes is further reduced because of the 
effect of perpendicular diffusion. However, in other simulations with different shock acceleration 
strength and parallel MEP (not shown here), if at the peak time of flux for W20 (W60), the flux 
for W20 (W60) is significantly different than that for E20, the reservoir phenomenon can not form 


the three observers. 


Reames et al. 


(1 1997 
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in normal diffusion coefficients. As a result, shock acceleration strength, parallel diffusion, and 
adiabatic cooling are the main factors in forming the reservoir phenomenon, and perpendicular 
diffusion is a secondary factor. 


3.3. Invariant Spectra Region 


There are some important characteristics in the invariant spectra region from our simulations 
(Figure [I])- If the observer is located at the eastern flank of the shock, the onset time of invariant 
spectra is earlier than the shock arrival. But if the observer is located near the central flank of 
shock, the spectra invariance begins approximately at the shock passages. Finally, if the observer 
is located at the western flank of shock, the onset time of invariant spectra are much later than the 
shock arrivals. From these results, we can better understand the invariant spectra region. 


Figure IS shows the invariant spectra region. In the picture, the gr een line is plotted by 


Reames et al. 


(| 19971) . and the red line is new in this work. According to 


Reames et al. 


19971) . the 


left side of the green line is the invariant spectra region, with the assumption that particles are 
quasi-trapped in the region behind the ICME, and the SEP fluxes gradually decrease as a result 
of parallel diffusion and adiabatic deceleration mechanisms, and in addition, there are also some 
leakage of energetic particles from ICME to eastern side of the upstream shock. In this sense, 
ICMEs play the pivotal role in the decay phase of fluxes. As a result, the invariant spectra region 
is determined by ICMEs’ propagation path plus some eastern side of upstream region. However, 
we suppose the invariant spectra region could be in the left side of the red line instead. In our 
simulations, ICME is not included, but in the propagation process perpendicular diffusion is 
included to reduce the spatial gradient in the fluxes, and adiabatic cooling is included to reduce 
the temporal variance. As the simulation results showed above, the spectra spatial and temporal 
invariance could result from the effects of shock acceleration strength, adiabatic cooling, and 
perpendicular diffusion. In this sense, it is possible that the invariant spectra region is not confined 
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by the ICMEs’ propagation path. Instead, the invariant speetra region eould be eonfined by the 
interplanetary shoek, but the region expands faster (slower) than the shoek at the eastern (western) 
flank, respeetively. 


4. DISCUSSION AND CONCLUSIONS 

We have studied interplanetary shoek aeeelerated SEPs propagation in three-dimensional 
IME The speetra observed by different observers are ealeulated, and the spatial and temporal 
invarianee in the speetra are reprodueed in the simulations. The following are our major findings. 

The adiabatie eooling effeet is the key faetor for forming temporal invarianee in the speetra. 
By ineluding the adiabatie eooling, for different energy ehannels, the flux deeay ratios are almost 
the same. The temporal invarianee results from the faet that all energy partieles deeay as the same 
ratios beeause of adiabatie eooling effeet. At the eastern flank of the shoek, the onset time of the 
speetra invarianee is earlier than the shoek arrival. Eor the eentral eases, however, the onset time 
of the speetra invarianee is elose to the time of shook arrival. At the western flank of the shoek, 
finally, the onset time of the speetra invarianee is later than shook arrival. In addition, the fluxes 
deeay muoh faster in the oases with adiabatie eooling. Without adiabatie eooling, the deeay phase 
of SEP fluxes are dominated by shoek aooeleration strength, parallel diffusion, and perpendioular 
diffusion, whioh are all varying with partieles’ energies. Therefore, the temporal invarianee does 
not exist without adiabatie eooling. 

Shook aooeleration strength, parallel diffusion, adiabatie eooling, and perpendioular diffusion 
are four important faotors in forming the spatial invarianee, whioh is in reservoir phenomenon in 
different energy ehannels. Shoek aooeleration strength parameters a and (f>c are set to 2 and 15° 
for 5 MeV protons in our simulations, respeetively. And these parameters also ohange with the 
momentum: a ~ p^'^, (pc ~ beeause the shoek aeeeleration strength deereases with higher 
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energy particles. Among the four factors, shock acceleration strength, parallel diffusion, and 
adiabatic cooling are the main factors in forming the reservoir phenomenon, and perpendicular 


diffusion is a secondary one. This conclusion is derived bas ed on our simu 


consistent with the observations. In 


Reames et al. 


(Il997h and 


Reames 


ations, and it is also 


(12013h . a gradual SEP event 


was detected by Helios 1, Helios 2, and IMP 8 during March 1, 1979 to March 11, 1979. In 
this event, the reservoir phenomenon appeared, and the onset time of SEPs observed by different 
spacecraft are very close because of the effect of perpendicular diffusion. At the peak time of flux 
observed by Helios 2 {IMP 8), the flux observed by Helios 2 {IMP 8) is close to that observed 


by Helios 1. The importance of the peak intensity of S EP flux o 


in forming the reservoir phenomenon is also noticed by 


Reames 


^ served by Helios 2 and IMP 8 


(12013h . however, the reservoir 


phenomenon is explained as a result of the disturbance of IMP caused ICMEs. In our model, 
the peak of flux is mainly determined by shock acceleration strength and parallel diffusion. 
Purthermore, the SEP fluxes decay as a similar ratio because of the effect of adiabatic cooling. At 
the same time, the latitudinal gradient in the SEP fluxes is further reduced because of the effect of 
perpendicular diffusion. Pinally, according to our model, the reservoir phenomenon appeared in 
this SEP event with the effects of shock acceleration strength, parallel diffusion, adiabatic cooling, 
and perpendicular diffusion. Observationally, shock acceleration s trength, diffusion coeffic ients. 


and adiabatic cooling change significantly in different SEP events (IKallenrode 


1996 


1997 ). As a 


result, the reservoir phenomenon can only form in some gradual SEP events with those controlling 
effect parameters in appropriate values. 

Based on our simula tions, a new inyarian t region is proposed. The new region is different 


from the one proposed by 


Reames et al. 


(119971) . There are two important characteristics in our new 


region. Pirst, if the observer is located at the eastern (western) flank of the shock, the onset time of 
temporal invariant in the spectra is earlier (later) than the shock arrival, respectively. Second, the 
spatial invariance in the spectra can also be formed without ICMEs. These two ch aracteristics ar e 


supported by observations, but are difficult to be explained in the previous model 


Reames et al. 
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(Il997h . 

In our model, we ignore the disturbanee of the IMF caused by ICME for the simplicity. In 
principle, the disturbance in the magnetic field can help particles redistribute in space. In future 
work, we intend to include a realistic three-dimensional ICME shock, so that the SEP acceleration 
and transport in the heliosphere can be investigated more precisely. 

The authors thank the anonymous referee for valuable comments. We are partly supported 
by grants NNSEC 41374177, NNSEC 41125016, and NNSEC 41304135, the CMA grant 
GYHY201106011, and the Specialized Research Eund for State Key Eaboratories of China. The 
computations were performed by Numerical Eorecast Modeling R&D and VR System of State 
Key Eaboratory of Space Weather and Special HPC work stand of Chinese Meridian Project. 
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Fig. 1.— Comparison of different energy protons observed by observers in the 1 AU eeliptie at 
different longitudes. The fluxes are normalized so their values are similar after all of them reaeh 
peaks. The vertieal lines indieate the shoeks passage of 1 AU. The adiabatie eooling is ineluded in 
simulations. 
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Fig. 2.— Same aslUexcept that the adiabatie eooling is not ineluded. 
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Fig. 3.— In the upper panel, eomparison of 5 MeV protons flux observed by the observers in 1 
AU eeliptie at different longitudes. The lower left and right panels show the speetra observed at 
different longitudes during time interval A and B, respeetively. The vertieal dashed lines indieate 
the shock passage of 1 AU. 
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The green line indieates the original speetra invariant region proposed by iReames et al. 
. The red line indieates a new speetra invariant region based on the simulation results in this 


paper. 













Table 1: Model Parameters Used in the Calculations. 


Parameter 

Physical meaning 

Value 

ysw 

solar wind speed 

400 km/s 

Vs 

shock speed 

870 km/s 

<Ps 

shock width 

60° 

a 

shock strength parameter 

2" 

4*c 

shock strength parameter 

15°'’ 

y 

injection spectrum 

5.5 

Be 

magnetic field strength at the Earth 

5 nT 

h 

particle radial mean free path 

0.2 AU" 

K± 

perpendicular diffusion coefRent 

0.1 xV 

no 

inner boundary 

0.05 AU 

ni 

outer boundary 

50 AU 


“for 5 MeV protons. 

*for 5 MeV protons. 

“for 5 MeV protons in the ecliptic at 1 AU. 
“^for 5 MeV protons in the ecliptic at 1 AU. 



